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ABSTRACT 



The Very Large Array (VLA) has been used in C-configuration to map an 
area ~ 0.3 deg 2 at 1.4 GHz with 5cr sensitivities of 0.305, 0.325, 0.380 and 0.450 
mJy beam -1 over four equal subareas. Radio properties are presented for 62 
detected sources. Deep optical imaging to Gunn r ~ 25 mag using the Hale 5-m 
telescope covering ~ 0.21 deg 2 is reported for a subset of 43 sources. This optical 
follow-up is much deeper than that of existing larger area radio surveys of similar 
radio sensitivity. Archival J, H and A-band photometry from the Two-Micron 
All Sky Survey (2MASS) is also presented. Using a robust likelihood-ratio 
technique, we optically identified 26 radio sources with probabilities ^ 80%, nine 
are uncertain/ambiguous detections, and eight are empty fields. Comparisons 
with a stellar synthesis model that includes radio emission and dust reddening 
suggests that the near-infrared-to-optical emission in a small, bright sub-sample 
is reddened by 'optically thin' dust with absorption Ay ~ 2 — 2.5 mag, 
regardless of morphological type. This is consistent with other, more direct 
determinations of absorption. The radio-to-optical(-near-infrared) flux ratios 
of early-type galaxies require significant contamination in the radio by an active 
galactic nucleus (AGN), consistent with the current paradigm. Using our 
simple modeling approach, we also discuss a potential diagnostic for selecting 
Ultraluminous Infrared Galaxies (ULIGS) to z ~ 1.6 from microJansky radio 
surveys. 



Subject headings: galaxies: active — galaxies: starburst — cosmology: 
observations — radio continuum: galaxies 
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1. INTRODUCTION 

The advent of deep radio surveys reaching flux densities well below 1 mJy (Mitchell 
& Condon 1985; Windhorst et al. 1985; Oort 1987; Windhorst et al. 1993; Hopkins et al. 
1998; Richards 2000) revealed a new population of faint sources more numerous than the 
AGN-powered radio galaxies dominating the strong-source population. This corresponds 
to a steepening of the differential source counts over non-evolving predictions at levels ^ 4 
mJy. The faint counts suggest that significant evolution has occured over the redshift range 
spanned by the observed population. Photometric and spectroscopic studies (Thuan et al. 
1984; Windhorst et al. 1985; Thuan & Condon 1987; Benn et al. 1993) suggest that the 
faint excess at 1.4 GHz is composed predominately of star- forming galaxies similar to the 
nearby starburst population dominating the Infrared Astronomical Satellite (IRAS) 60-/xm 
counts (Benn et al. 1993). Indeed, this is supported by the strong correlation between 
radio (1.4 GHz) and far-infrared (60-yU.m) flux densities of disk galaxies (Helou, Soifer & 
Rowan- Robinson 1985), implying a significant proportion of starburst galaxies at faint radio 
flux densities. 

The overall observed source-count distribution from faint (/xJy) to bright (Si. 4 ^ lOmJy) 
flux densities cannot be explained by starbursts alone. Evolutionary models of radio source 
counts need to invoke two separate populations (eg. Danese et al. 1987; Rowan-Robinson 
et al. 1993; Hopkins et al. 1998). Condon (1989) describes these populations as starbursts 
and monsters, each powered by different mechanisms: 'starbursts' deriving their radio 
emission from supernova remnants and HII regions and 'monsters' from compact nuclear 
related activity (eg. active galactic nuclei; AGNs). The proportion of AGN is much greater 
at higher flux densities S\ A ^ 10 mJy (Kron, Koo & Windhorst 1985; Gruppioni et al. 
1998), where a majority are associated with classical radio galaxies exhibiting extended 
(FRI and FRII-type) morphologies (Fanaroff & Riley 1974). The optical counterparts to 
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sources at bright radio flux densities 51.4 ^ 1 m Jy is comprised mostly of ellipticals while 
at sub-mJy to /iJy levels, the optical counterparts are identified as blue galaxies exhibiting 
peculiar (compact, interacting and irregular) morphologies (Kron, Koo & Windhorst 1985; 
Gruppioni et al. 1998). Studies of faint radio sources, namely their stellar population, 
how they evolve with redshift, and how they relate to local normal galaxies is progressing 
rapidly, however, much remains to be learned from the faintest (/iJy) radio populations at 
redshifts 1 ^ z <^ 2. 

Radio surveys are insensitive to the effects of absorption by dust which is known to bias 
surveys severely at optical/UV wavelengths. This is particularly important for derivations 
of the cosmic history of star formation and its relation to hierarchical models of galaxy 
formation. Optical/UV studies have shown that there is an increase in both the space 
density of star-forming, morphologically disturbed galaxies (eg. Richards et al. 1999) and 
also, the global star formation rate with redshift to z ^ 1 (eg. Madau et al. 1996). Similar 
evidence is emerging from studies of sub-mm sources (Blain et al. 1999) and amongst the 
faint radio population at 1 <^ z <^ 2 - the redshift range probed by the deepest surveys (eg. 
Cram 1998; Haarsma et al. 1999). 

There has been much speculation as to whether global star-formation rates (SFRs) 
derived from radio observations exceed those determined from optical/UV studies. Cram et 
al. (1998) note that systematic discrepancies may exist between the various star formation 
indicators, which are not well understood (see also Schaerer 1999). It is encouraging to see 
however that Haarsma et al. (2000) find global SFRs derived from 1.4GHz observations of 
the Hubble Deep Field to exceed optically detetermined values by a factor of a few out to 
z ~ 1. Indeed, an analysis of Balmer decrements and optical-near-infrared colors in star 
forming galaxies by Georgakakis et al. (1999) from the Phoenix Deep Survey (Hopkins et 
al. 1998), finds evidence for visual extinctions from one to a few magnitudes. Currently, 
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about 20% of existing micro-jansky radio samples remain unidentified to / = 25 mag in 
Hubble Deep Field images (eg. Richards et al. f999). A majority of these could represent a 
significant population of dust-enshrouded starbursts and/or AGN at high redshift. These 
results are in support of efforts to further understand the dust properties of star-forming 
galaxies. 

The primary aim of this paper is to introduce a new complete sample of radio sources 
selected at 1.4 GHz, uniformly selected over the flux range S ^ 0.3 mJy (5a) from an area 
covering ~ 0.3 deg 2 . Although much larger area surveys to deeper radio sensitivities have 
been carried out (eg. Hopkins et al. 1998), the present study reports the results of more 
sensitive optical observations. Archival near-infrared data for a subset of the sample are 
also presented. The near-infrared data is from the ongoing 2MASS project, and represents 
a unique aspect of this study in the identification of radio-selected starbursts. Although we 
currently lack valuable spectroscopic information, we combine radio-near-infrared-optical 
flux ratios, radio maps, and optical images to explore the properties of the entire sample. 
Our deep optical identifications provide the opportunity to asses the importance of dust 
in star-forming systems via the observed radio-to-optical and near-infrared-to-optical 
colors. Simple stellar synthesis models that include radio emission and reddening are used 
to constrain possible amounts of absorption. 

This paper is organised as follows. In §2, we discuss the radio observations and data 
reduction, present the radio catalog, and compare our results with data available from 
(shallower) all sky radio surveys. The optical photometric observations, their reduction and 
astrometric calibrations are discussed in §3. Our method for radio-optical identification, 
the archival near- infrared data, radio-optical image overlays, and our optical/near-infrared 
catalog are presented in §4. A study of the radio, near-infrared, and optical colors and 
constraints on synthesis models incorporating dust is presented in §5. An application of our 
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color-color analysis to select high-redshift ultraluminous infrared galaxies from deep radio 
surveys is discussed in §6. All results and future prospects are summarised in §7. 



2. RADIO OBSERVATIONS 
2.1. Strategy 

Observations were made with the VLA C-configuration at 1.4 GHz on 1998 December 
19. This configuration yields a good compromise for resolution and surface-brigntness 
sensitivity. The 5a confusion limit for this configuration is only 50/iJy/beam at 1.4 GHz 
(Mitchell Sz Condon 1985) since the synthesised beam size (Full Width at Half Power; 
FWHP) is ~ 15 arcsec. The resulting radio positions have rms errors ~ 1 arcsec, except for 
extended sources with multiple components, sufficient for making optical identifications. At 
1.4 GHz the FWHP of the VLA primary beam is 31 arcmin and approximately corresponds 
to the diameter of our final imaged field. This relatively large coverage avoids field-to-field 
variations in source counts induced by high-redshift clustering. Although surveys at higher 
frequencies (eg. 8 GHz) can reach lower flux densities than at 1.4 GHz, most radio sources 
have spectral indices a ~ 1 (S oc v~ a ) so that the population being sampled is similar. 

To provide uniform sensitivity over the full area of our field, we observed seven positions 
arranged in a filled hexagonal pattern with a separation of 26' between pointing centers. 
The resulting root mean square (rms) map noise is thus nearly constant (cf. Condon et al. 
1998). Our field was centered on RA(2000)=00 /l 13 m 12 s .O, Dec.(2000)=+25°54'44". This 
field was chosen for its relatively low foreground galactic-cirrus emission likely to affect 
optical/near-infrared identifications and also for the absence of bright radio galaxies. The 
integration time on each pointing was ~ 1 hour. This allowed us to reach an rms noise of 



60yuJy in regions free from bright contaminating sources (see § [2.3| for more details). 
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Our observations were made in spectral line mode with four Intermediate Frequencies 
(IFs), each divided into 7 spectral channels of width 3.125 MHz. The advantage of this 
mode is to minimize bandwidth smearing (i.e. chromatic aberration) which reduces the 
point-source sensitivity away from the pointing center and causes appreciable image 
distortion. Additionally the spectral line mode is less prone to narrowband interference 
noise spikes. With continuum mode however, we would have had a little over twice the 
bandwidth and a factor a/2 lower noise. 

2.2. Data Reduction 

The data were analysed with the NRAO MPS reduction package. We observed 3C 
48 to calibrate the visibility amplitudes, using S = 16.5 and S = 15.9 Jy at 1.365 and 
1.435 GHz, respectively. The calibration was applied using the SPLIT task. The calibrated 
data from each pointing were edited and imaged separately, CLEANed, and restored with 15 
arcsec FWHP Gaussian beams. The seven separate images were weighted and combined 
as described in Condon et al. (1998) to produce a final 33' x 33' map with nearly uniform 
sensitivity and corrected for primary-beam attenuation. 

2.3. Noise and Source Extractions 

The resulting rms noise of our final map after correcting for primary-beam attenuation 
is not uniform over the entire field, but increases by up to 30% near a strong (^ 400 mJy) 
source near the edge of our field. Despite this variation in sensitivity, we were able to divide 
our 33 x 33arcmin 2 field into four equal (~ 16.5 x 16.5arcmin 2 ) regions within which the 
rms noise varies by no more than a few percent. These constant- noise regions simplify the 
application of an automated source-finding algorithm over a single continuous region (see 
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below). The lowest and highest rms noise amongst these regions was ~ 60.9 yuJy/beam and 
~ 90.3 //Jy/beam respectively. See Table [3] for the region definitions. The rms noise of each 
region was estimated from the Gaussian core of the amplitude distribution of the pixel 
values as produced by the AIPS task IMEAN. In Figure p], we show the distribution in pixel 
values of our entire 33' field. The rms deviation in peak flux density derived from a fit to 
this histogram is ~ 69/iJy/beam. Figure ^ shows a contour map of our entire radio field. 

Each constant-noise region in Table [I] was used for the source extractions. Within each 
region, we searched for radio sources down to a peak flux density > 5 times the rms value 
of the region. The sources were extracted by the AIPS task SAD which uses Gaussian fits 
to estimate the fluxes, positions, and angular sizes of the selected sources. For faint sources 
however, unconstrained Gaussian fits may be unreliable (see Condon 1997). For this reason, 
we adopted the following method for the source extraction: first, we ran the task SAD with 
a 3<r rms detection threshold to obtain an initial list of candidates, we then derived the peak 
flux densities of the faint sources (with 3a < S pea k < 7a) using the MAXFIT task. This task 
uses a second order interpolation algorithm and is known to be more accurate. Only sources 
with a MAXFIT peak flux density > 5a rms were retained. For these faint sources, the total 
flux density was estimated using the task IMEAN, which integrates the (median background 
subtracted) pixel values in a specific rectangle. The rectangle was chosen to encompass 
as much of the source as possible. For all other parameters (sizes, positions and position 
angles) we retained the values obtained from the initial Gaussian fits. Only two sources 
had irregular morphologies showing multiple components. For these, the total (background 
subtracted) flux was determined using the task TVSTAT which allows an integration over a 
specific irregular area defined to encompass as much of the source as possible. 

The numbers of sources detected in each constant-noise region are summarised in 
Table jl} A total of 62 sources were detected with flux densities > 5a rms over an area of 
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Fig. 1. — Logarithmic histogram of pixel flux densities in our complete radio map. 
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Fig. 2. — Contour plot of our radio map covering an area 33 x 33 arcmin 2 = 495 x 495 pixels. 
Contour levels are 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6 mJy/beam. Note the increase in rms 
noise from the bright off-field source in the top right corner. 



- 11 - 



0.303 deg 2 . Within poisson uncertainties, this number is consistent with source counts 
from surveys by Ciliegi et al. (1999) and Hopkins et al. (1998). To our limiting (mean) 
sensitivity of 0.35mJy and within 0.303 deg 2 , they report a source count of typically 70 ± 8. 
This confirms the reliability of our radio source detections and flux density estimates. 

Table 2 shows the full radio catalog which reports (in column order): the source name; 
RA and Dec.(J2000); errors in RA and Dec; the peak flux density Sp; error in Sp; the 
total flux density St] error in St] the full width half maximum (FWHM) of the major 
and minor axes 9 m arid 9 m (determined from Gaussian fits); the position angle PA of the 
major axis (in degrees); and the rms errors associated with 9m, 9 m , and PA respectively. 
The different components of multiple sources are labelled "CI" and "C2" . In Figure we 
show the distribution of peak flux densities and the total to peak flux ratio as a function 
of peak flux for all sources. Sources with ratios St/Sp < 1 in Figure |3|b are primarily 
due to uncertainties on measured fluxes as estimated from the two dimensional Gaussian 
fits. Uncertainties in peak fluxes are typically <^ 10%, while total flux estimates are more 
uncertain due to a sensitive dependence on the Gaussian fitting procedure. Errors in total 
fluxes are typically ^ 30%. Contour maps of radio sources with available optical data are 
shown in Figure 

2.4. Comparison with the NVSS Catalog 

The 33 x 33 arcmin 2 region that we observed with the VLA was also covered by the 
NRAO VLA Sky Survey (Condon et al. 1998, NVSS). The NVSS covers the sky north of 
5 = —40° at 1.4 GHz with 45" resolution and limiting flux density of ~ 2.25 mJy (5a). To 
this limit, 17 of our sources were found to be in common with the NVSS public catalog. 
One source in our catalog however is a double component source and is unresolved by the 
NVSS. A comparison of flux densities derived in this study with those from this catalog 
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Fig. 3. — (a) Distribution of peak flux density and (b) total to peak flux ratio as a function 
of peak flux. 



-13- 

is shown in Figure It is evident that our derived flux densities are on average slightly 
lower than those from the NVSS. This was also reported in a larger comparison study 
by Ciliegi et al. (1999) using a similar observational set-up and can be explained by the 
difference in resolutions used in the two surveys: 15" vs. 45". High-resolution surveys tend 
to miss flux from low-surface-brightness emission. Although the effect is only marginal for 
bright sources, it may become important for attempts to detect faint, low-surface-brightness 
objects at high redshift. 

3. OPTICAL PHOTOMETRIC OBSERVATIONS 

3.1. Observations and Data Reduction 

Optical CCD photometry of our radio field was carried out on the 5-m Hale telescope 
at Palomaif] Observatory during the nights of 27-29 August 1997. These fields were initially 
selected for subsequent deep mid-IR imaging with the Wide Field Infrared Explorer mission 
(WIRE), but the mission failed to perform to expectations. The Carnegie Observatories 
Spectroscopic Multislit and Imaging Camera (COSMIC; Kells et al. 1998) mounted at the 
f/3.5 prime focus with a TEK 2K CCD was used to image nine 9.7 x 9.7arcmin 2 fields in 
the Gunn r (6550A) filter. Each optical field comprised of three optical pointings offset 
by 2', each with integrations of 600 sec giving a total 1800 sec per field. This resulted in 
a limiting magnitude of r ~ 25 mag (5a). The seeing was typically 1 — 1.4" (FWHM). 
The optical fields do not cover our entire 33 x 33 arcmin 2 radio field. The nine slightly 
overlapping optical fields correspond to an areal coverage 27.5 x 27.5 arcmin 2 , or about 70% 
of the radio map. 



5 Operated by California Institute of Technology, Pasadena, CA 91125 
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Fig. 4. — The peak flux density obtained with our survey versus that reported for the same 
sources in the NVSS public database. 
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The CCD data were reduced with standard tasks in the IRAF package. Frames were first 
bias-subtracted, then flat-fielded using dome flats. Bad pixels and columns were removed 
by interpolating between adjacent pixels and lastly, the individual dithered frames were 
median combined to remove cosmic ray hits. Calibration was performed using standard 
stars from Thuan & Gunn (1976). These were used to correct for atmospheric extinction 
from varying airmass and provide the instrumental zero point. Photometric uncertainties, 
estimated using these standards, are no more than ps 0.05 mag. 

Sources in the reduced optical frames were extracted using the DA0PH0T package in IRAF 
(Stetson 1987). This package has the benefit of performing photometry in crowded fields, 
which is the case in most regions of our optical fields. It performed the following steps: 
first, sources were extracted above a given threshold, given as a multiple of the total CCD 
noise (sky and read noise, <J tot ). We adopted threshold values of 4.5-5a tot , just high enough 
to avoid large numbers of spurious detections. Second, simple aperture photometry was 
performed on these identified sources. This required a specification of the aperture radius 
which is likely to contain most of the light of our target source, and width of a surrounding 
annulus to estimate and subtract the sky background. We adopted a radius of six arcsec 
and annulus width of four arcsec. Next, a point spread function (PSF) was determined in 
each of our nine fields. This involved an interative technique to remove contamination from 
neighbouring sources in crowded fields near our PSF candidates. Simultaneous PSF fitting 
on all sources was then performed to identify sources which were previously hidden in 
crowded regions. Finally, the magnitudes determined from PSF fits were aperture corrected 
to a common aperture size as used on our standard stars. Aperture corrections were 
typically 0.22 mag. A final visual inspection removed any spurious detections. A total of 
w 300-390 sources were extracted from each 9.7 x 9.7arcmin 2 field to a limit of r ~ 25 mag. 
Previous optical surveys find typically 380-520 sources within this area to this limit, and 
the variation is primarily due to clustering. Such fluctuations are found to be significant on 
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such relatively small scales (eg. Metcalfe et al. 1991). 



3.2. Astrometry 

The astrometry on the optical images was based on 10-12 APM catlog stars in each 
field (Maddox et al. 1990). The ccmap and cctran tasks in the IRAF immcoords package 
were used to compute plate solutions relating pixel positions to astrometric coordinates. 
Astrometric coordinates for all sources on the frames were then determined. By comparing 
the positions of several sources common to the APM catalog and our fields, we found that 
our rms position uncertanties are typically ^ 0.9". 



4. OPTICAL /NEAR-INFRARED IDENTIFICATION OF RADIO 

SOURCES 

4.1. Method for Optical Identification 

We assigned optical identifications and estimate their reliability using a robust 
likelihood ratio (LR) analysis. This general method has frequently been used to assess 
identification probabilities for radio and infrared sources (eg. de Ruiter, Willis & Arp 
1977; Prestage & Peacock 1983; Sutherland & Saunders 1992; Lonsdale et al. 1998). The 
method, which computes the probability that a suggested identification is the 'true' optical 
counterpart, is outlined as follows: 

For each optical candidate % in the search area of some radio source j, we calculated 
the value of the dimensionless difference in radio and optical positions: 



r, ij 



(at - OLjf (5i - 5j) 2 



1/2 

(1) 



where the cc's and <5's represent right ascensions and declinations respectively, and cr's 
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standard deviations. We chose a moderately large search radius of 10" to allow for the 
maximal position uncertainties: a opt ~ 1.4" and <r ra d ~ 1.5" (assuming 5(T e ff, where 
a e ff = (<7o pt + o"rad) 1//2 )- Such a radius is also small enough to avoid large numbers of chance 
associations. 

Given r^, we must now distinguish between two mutually exclusive possibilities: (1) 
the candidate is a confusing background object that happens to lie at distance from 
the radio source or (2) the candidate is the 'true' identification that appears at distance 
Tij owing solely to radio and optical position uncertainties. We assume that the radio and 
optical positions would coincide if these uncertainties were zero. This assumption however 
is not valid when the centroid of an extended radio source is used, and is further discussed 
below. 

To distinguish between these cases, we compute the likelihood ratio LR^, defined as: 

hfi/2 



exp 

LR^ = ■ — — -, 1/2 - 



(2) 



where N(< raj) is the local surface density of objects brighter than candidate i. The 
likelihood ratio LRij is simply the ratio of the probability of an identification (the Rayleigh 
distribution: r exp (— r 2 /2)), to that of a chance association at r (27riV(< mi)a a as). LR^ 
therefore represents a 'relative weight' for each match ij, and our aim is to find an optimum 
cutoff LR C above which a source is taken to be a reliable and likely candidate. Its advantage 
over alternative methods (purely based on finding the lowest random chance match, eg. 
Downes et al. 1986) is that it allows for a possible distant candidate to still be the 'true' 
identification even when there is still a high chance of it being a spurious background 
source. 

It is important to note that our form for LR^ (eq. 0) slightly differs from that used 
by earlier studies (eg. Lonsdale et al. 1998) in that it doesn't contain the multiplicative 
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"Q" factor in the numerator. This factor represents the apriori probability that a "true" 
optical counterpart brighter than the flux limit exists amongst the identifications. For our 
purposes, we will treat LRij as simply a a relative weight measure for each radio-optical 
match, just for the purposes of assigning an optimal cutoff for reliable identification (see 
below). We are not concerned with its absolute value, which is required when computing 
formal probability measures from LR. For simplicity, we have therefore set Q = 1 in this 
work. 

The optical surface density as a function of magnitude to be used in computing LR was 
determined from the total number of objects visible in our optical frames. The variation in 
surface density in the vicinity of each radio source caused by possible clustering effects was 
found to be small: no more than 5% on 2 — 3' scales. 

The distribution of LR values for all possible radio source-candidate matches is shown 
by the shaded histogram in Figure |5]. Following Lonsdale et al. (1998), we generate a truly 
random background population with respect to the radio sources by offseting the radio 
source positions by « 30". LR values for each radio source were then re-computed and their 
distribution is given by the thick-lined histogram in Figure [5]. A comparison of the number 
of associations for (true) radio source positions with the number of associations found 
for random (offset) positions will enable us to determine a critical value LR C for reliable 
identification. From these distributions, we compute the reliability as a function of LR: 



where N trU e and N ran d om are the number of true and random associations respectively. The 
reliability computed in this way also represents an approximate measure of the identification 
probability for a candidate with given LR. 



reliabilities are ^ 85% because few random associations exceed this value of LR (Fig. [|). As 




(3) 



Figure |G] illustrates the reliability as a function of LR. Above log (LR) ~ 0.5, the 
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Fig. 5. — Distribution of LR for radio-optical matches at "true" radio positions (shaded 
histogram), and at "random" radio positions (thick-lined histogram). 
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Fig. 6. — Reliability as a function of LR. See equation (§). 



-21 - 



a good working measure we therefore assume a cutoff log (LR C ) ~ 0.5 above which a source 
is taken to be a likely candidate. The determination of reliabilities via the LR method 
is insensitive to variations in N(< mi) across the field or uncertainties in its derivation, 
and also the assumption of Gaussian error ellipses in the radio and optical positions. Such 
uncertainties are "normalised out" when one computes the ratio of random to true number 
of associations within a search radius when estimating the reliability (eq. ^]). Lonsdale et 
al. (1998) have shown that the absolute value of LR itself depends on the characteristics 
of the source population being identified (eg. stars versus galaxies). Different populations 
(assuming they could be classified a-priori using some diagnostic) map into different 
underlying surface densities at the 'identifying' wavelength, implying that distributions 
in LR (eq. 0) will also be different. For a robust determination of the reliability in such 
situations, see Lonsdale et al. (1998). 

There are two complications to consider in the above method. The first is when one 
attempts to identify extended (or resolved) radio sources with this method. For all radio 
sources, we have used the positions of centroids derived from two-dimensional Gaussian 
fits in computing the LR for optical candidates. For unresolved sources with (Gaussian 
fitted) sizes <^ 15" (the synthesised FWHP beam width), the source is likely to have a 
compact central component and the optical position is expected to lie close to its quoted 
radio centroid. For an extended (resolved) source however, the radio and optical positions 
may differ considerably since errors in the radio centroids are only <; 2". In such cases, if 
LR < LR C the identification may still be valid, since its low LR value could purely be due 
to a real large positional offset. The second complication is when a radio source has more 
than one optical candidate within its search radius with LR > LR C . This occured in about 
20% of cases and was primarily due to contaminating stars. We assess these ambiguities and 
increase the robustness of our identifications by visually examining all optical candidates 
according to the following criteria: 
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1. If candidates have Li? < Li? c for a radio source with 9 min ^ max ^ 15", then identification 
is classified as uncertain. 

2. If candidates have very low reliability, LR <C LR C (for unresolved radio sources), or 
there are no objects in the search radius, then radio source is classified as empty field. 

3. If LR <^ LR C , i.e. where reliability is moderately "low", then identification is also 
uncertain. 

4. If more than one optical candidate exists with LR > LR C , then only source(s) with 
extended (galaxy-like) optical profile is taken as the identification. Point sources 
associated with quasar nuclei are not considered in our identification scheme due to 
their relatively low surface density compared to galaxies (~ 1 : 4000) in sub-millijansky 
radio samples. 

5. For unique, LR > LR C candidates, its optical profile is also checked for confirmation. 

4.2. Results 

Of the 62 radio sources, 43 lie within our optically imaged 27.5 x 27.5arcmin 2 field. 
We found optical identifications for 26 to r ~ 25 mag with reliabilities R i( i ^ 80%. Four 
sources have identifications classified as uncertain owing to a moderately low identification 
reliability of Rid <^ 78% (and log (Li?) <^ 0.4). Five more are uncertain because they 
have extended radio structure and large possible positional offsets between optical and 
radio centroid positions. Eight radio sources lie in "definite" optical empty fields with no 
candidates brighter than r ~ 25 mag. 

Other optical follow-up studies found similar results. Georgakakis et al. (1999) 
identified as 47% of sources to R — 22.5 mag from the Phoenix Deep Survey (Si.4gh z > 0.2 
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mJy) (Hopkins et al. 1998). Deeper indentificatioris of sources as faint as Si.4GHz — 40^ Jy 
from Hubble Deep Field images revealed a 80% success rate to / = 25 mag (Richards et al. 
1999). Ignoring the uncertain identifications in our study (from criteria 1 and 3 above and 
which are excluded from our analysis), we find that «18% of our sources are unidentified to 
r ~ 25. Accounting for differences in bandpasses and sensitivity, this is broadly consistent 
with the above studies. Figure |7| shows the distribution of apparent magnitude r for all 
reliable (robust) and uncertain identifications in our sample. 



4.3. Near-Infrared Data 

Near-infrared data in the J(1.25/xm), if(1.65/iin) and K s (2.17fim) photometric 
bandpasses were obtained from the Two Micron All Sky Survey (2MASS) project database. 
For multi-band detection of point sources, this survey is currently scanning the sky to 
sensitivities 16.5, 16.0 and 15.5 mag at signal-to-noise ratios ~ 7, ~ 5 and w 7 respectively 
in J, H and K s . The data relevant to this study are not yet released in the public catalogs, 
and was retrieved from the 'internal working database' at IPAC^ Photometry in this 
database was determined using custom PSF-fitting software and algorithms are described 
in Cutri et al. (2000). Since such data have not been subjected to the rigorous quality 
assurance as that in the public release catalogs, we have examined individual images for 
quality and any possible systematic uncertainties in the photometry. 

To maximise the possible number of detections, we searched for near-infrared 
counterparts to each radio source with a conservatively low signal-to-noise ratio threshold 
of ~ 5 in each band. In cases where a source was detected in only one or two of the 
three bandpasses: J, H and K s , we note its "band- filled" 95% confidence upper-limit in 



6 The Infrared Processing and Analysis Center, California Institute of Technology. 
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Fig. 7. — The distribution of apparent r-band magnitudes for the 26 identified radio sources 
(open histogram) and nine uncertain identifications (filled histogram). Eight empty fields 
are represented by the r > 25 magnitude bin with arrows. 
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the undetected band. In other words, the 2MASS catalog also reports upper limits to 
measurements in an undetected band by placing an aperture over the position inferred from 
detections in other bands. The image pixel scale of the 2MASS detectors is 2.0" and the 
positional uncertainties are <^ 0.5". Owing to the relatively shallow flux limits of the 2MASS 
survey, the background source surface density is low enough that chance associations with 
radio positions are very unlikely. 

We searched the 2MASS database for near-infrared counterparts to our 43 radio sources 
that have available optical information from our deep optical survey and examined their 
images for quality. We found 7 reliable matches with 6 detections in J, 7 in H, and 4 in K s 
at the ^ 5<T level. Upper limits were available for the remainding "band-filled" values. 

Results of our optical identification analysis and available near-infrared data are shown 
in Table 3. In column order, this table reports: the radio source name (see Table 2); 
RA and Dec.(J2000) of the optical counterpart of the radio source; optical-radio position 
separation (S ra( i-o P t) i n arcsec; logarithm of the liklihood ratio {LR); reliability of the optical 
identification (see eq.||); apparent r-band magnitude; J, H and K s magnitudes with errors 
or 2<r upper limits; r — H color; r — K s color, and last, the optical morphology if the optical 
counterpart is visually resolved with size ^ 59fwhm(PSF). 

Optical morphologies were determined from light profile fitting of individual sources 
and comparison with i? 1//4 -law (elliptical-like) and exponential (disk-like) profiles. In cases 
where a disturbed or interacting morphology is apparent, then it is designated to have an 
irregular (labelled as 'Irr') morphology. Sources unresolved in the optical with typically 
59fwhm(PSF) are designated as 'unknown' and labelled as '?' in Table 2. 
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Fig. 8. — Optical image (r-band) and radio contour overlays for the 43 sources with available 
optical information. Contour levels are la, 2a, 3a, 4a, 5a, 10a, 20a, 30a, 40a, 50a, lOOcr and 
200a, where a refers to the local rms noise (see Table H). Optical candidates are indicated 
within broken lines — see Table 3 for reliability estimates. Maps labelled with "E" in upper 
left corner represent optical empty fields. 
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4.4. Optical and Radio Map Overlays 

In Figure § we show the optical image - radio map overlays for the 43 radio sources 
with available optical information. A visual inspection of the optical images of resolved 
counterparts shows a diverse morphological mix, consistent with previous studies. About 
40% of our optical identifications have elliptical/disk-like morphology, while ~ 10% can 
be identified as exhibiting peculiar (either interacting or disturbed) morphologies. It is 
important to note that these 'disturbed' sources are based on visual inspection alone 
and their morphology could still be uncertain until future spectroscopy confirms their 
nature. The elliptical/disk hosts also tend to be associated on average with sources of 
relatively brighter radio flux density 2 mJy) than the irregular class. This is consistent 
with previous radio-optical identification studies which find an increasing fraction of 
irregular-type galaxies at S1.4 ^ 2 mJy (eg. Kron, Koo, & Windhorst 1985, Hammer et 
al. 1995, Gruppioni et al. 1998) and a significant number of elliptical galaxies hosting the 
brighter extended radio galaxies and AGN (Condon 1989). 

A further observation is the unique radio structure exhibited by our eight optical 
empty field sources with r > 25 mag. These are represented by maps in Fig. || labelled 
by the letter "E". All show compact and symmetric (presumably unresolved) structures 
and could represent either of the following: distant (possibly dusty) AGN where with our 
radio sensitivity, we could have detected a nominal FR-I galaxy to z ~ 1.3, or, nearby, 
compact dusty starbursts at z <^ 0.3 as constrained by typical starburst luminosities: 
Liaghz ^ 10 23 WHz _1 . The second explanation for the nature of the empty fields is more 
plausible, given that the majority of submillijansky radio sources are associated with 
star-forming galaxies and less than 5% are usually identified with bright FR-Is at z ^ 1 
(Kron, Koo & Windhorst 1985). 

EDITOR: PLACE FIGURE | HERE. 
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5. ANALYSIS OF RADIO AND OPTICAL-NEAR-INFRARED COLORS 

This section presents an analysis of flux ratios between the radio, near-infrared, 
and optical bands to explore possible contributions from AGN and starbursts to the 
observed radio emission as well as the importance of absorption by dust. Because we lack 
spectroscopic information, our analysis treats the sub-mJy sources as one homogeneous 
population and uses a simple stellar synthesis model to interpret its properties quantitatively. 

5.1. A Simple Synthesis Model 

We can predict the radio-to-optical (-near-IR) flux ratios and r — K colors for a range 
of galaxy types using the stellar population synthesis code of Bruzual & Chariot (1993) 
(hereafter BC). On its own however, the BC model does not directly predict the amount of 
radio emission expected from a star-forming galaxy, which we need for the determination of 
flux ratios involving the radio band. We do this by relating the star-formation rates derived 
from empirical calibrations involving the UV and radio bands as follows. 

The 1.4 GHz radio emission from star-forming systems is believed to be primarily 
synchrotron emission from cosmic rays accelerated in supernova remnants plus a small 
(~ 10%) thermal contribution from HII regions (Condon & Yin 1990, Condon 1992). Thus 
to a good approximation, the radio luminosity is taken to be proportional to the formation 
rate of stars with M > 5Mq: 

SFR(M > 5M W ) = f; L4GHz i-Afoyr- 1 . (4) 

V ~ UJ 4 x 10 28 erg s^Hz -1 U V ; 

(Condon 1992). These same massive stars will also contribute significantly to the UV 

continuum emission in the range ~ 1200-2500 A. In particular, there have been many 

different calibrations of the SFR from the UV-flux. For a Salpeter initial mass function 

(IMF) from ttlMq to IOOMq, the calibration of Madau et al. (1998) (which assumes no 
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dust correction) yields 

sfr ( m > mM@) = Qm ( „ 4xl0 t v gs - 1Ifa -- ) M ° yr " (5 > 

We have modified the initial relation of Madau et al. (1998) to include the factor Q m , which 
represents the fraction of stellar masses contributing to the observed SFR, 

_ Lm® Mip(M)dM 

Qm — 10OM0 ' (^) 

f 0AM O Mip(M)dM 

where ip(M) oc M _x is the IMF. For m = 0.1, we have Q m = 1- Assuming a Salpeter IMF 
(x = 2.35), we find that for stars with M > 5Mq, Q m ~ 0.18. With this fraction, and 
equating the two SFR calibrations (eqs 0] and 0), we find that the luminosity densities at 
1.4 GHz (L 1A chz) and ~ 2100A (L uv ) are very nearly equal. We therefore assume that 
the rest frame 1.4 GHz flux density is directly given by the flux density at ~ 2100A as 
specified by the synthesized model spectrum. In general terms, the observed radio flux (in 
the same units as the synthesised UV spectrum) can be written: 

f u (lAGEz) obs = (1 + z) 1 - Q / l/ (2100A) rest , (7) 

where a ~ 0.8 is the radio spectral index (Condon 1992) and f l/ (2100A), rest is the rest frame 
(unreddened) UV spectral flux. We must emphasise that this relative radio flux is only 
that associated with the star-formation process. Possible additional contributions, such as 
contaminating AGN, are not considered in this model. 

We calculated flux ratios involving the radio, near-infrared and optical bands using 
evolutionary synthesis models for ellipticals (E/SO), early (Sab/Sbc) and late (Scd/Sdm) 
type spirals, and "very blue" starbursts (SB). These are meant to represent the possible 
contributions to the sub-mJy radio sources, and each class is defined by a characteristic 
star formation rate as a function of time. As supported by local observations (eg. Gavazzi 
& Scodeggio 1996), we assumed that E/SO and Sab/Sbc galaxies have an exponentially 
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decaying SFR of the form ip(t) oc r _1 exp (— t/r), where r is the e-folding time. Values of 
t = 1 and r = 8 Gyr were adopted for the E/SO and Sab/Sbc galaxies respectively For 
late-type spirals (Scd/Sdm) and young starbursts (SB), we assumed constant SFRs with 
different ages. All models used to generate the spectral energy distributions (SEDs) are 
summarised in Table f|. The models assume H = 50kms _1 Mpc -1 , go = 0.5 and a galaxy 
formation redshift Zf = 10, which corresponds to an age of 12.7 Gyr at z = 0. 

To explore the effects of dust on our flux ratios and colors, each model SED was 
reddened in the source rest frame with an extinction curve £(A) = A\j Ay characteristic 
of the Small Magellanic Cloud (SMC), which appears to be a good approximation for the 
ISM of nearby galaxies (Calzetti et al. 1994). This approximation is most accurate for 
the reddest wavelengths of starburst galaxy spectra (7000A-3/im), although breaks down 
severely at A ^ 2500A(Calzetti et al. 1999). We have used the analytical fit for £(A) as 
derived by Pei (1992) for the SMC. For simplicity, we assumed that the dust is distributed 
in a homogeneous foreground screen at the source redshift. 

5.2. Data and Model Comparions 

In Figure [| we plot the radio-to-optical flux ratio, R(lA/r), defined as 

R(lA/r) = log (5 L4 /mJy) + 0.4 r, (8) 

where Si. 4 and r are the radio flux and optical r-band magnitude respectively. The 
distribution seen in observed values of R(lA/r) is larger at the faintest optical magnitudes 
r > 21.5. There are few galaxies however with r < 21.5, and nonetheless, the scatter at 
r > 21.5 is consistent with that found at r ^ 19 in a larger follow-up study of sub-mJy 
radio sources by Georgakakis et al. (1999). 

Figure [U] shows R(lA/r) as a function of r — H color for all sources with available 
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Fig. 9. — Radio-to-optical (r-band) flux ratio (see eq. [||) as a function of r magnitude for 
all radio sources with available optical information. For the observed range of r-magnitudes, 
the solid line represents the prediction at the limiting flux of the survey Si.4gh z ~ 0-3 mJy. 
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optical and near-infrared data. Our reason for using r — H color is that the H band 
yielded more "definite" detections than the other near-infrared bandpasses. Although the 
numbers are still relatively small, this facilitates the best comparison with the synthesis 
models. The predictions for four galaxy types (see Table |] and §5.1) are shown for no dust 
reddening (thin curves) and a rest-frame extinction Ay = 2.5 mag. The morphological 
mix of data shows a relatively large dispersion in r — H color that is more consistent with 
the range predicted by the models that include dust. This suggests that on average, the 



optical-to-near-infrared continua of most sources in Fig. [Il] are reddened by a uniform 
(possibly "optically thin") dust component with Ay ~ 2 — 2.5 mag absorption. This 
measure is consistent with spectroscopic studies of nearby starbursts by Calzetti et al. 
(1996), Meurer et al. (1997) and photometric modelling by Nakata et al. (1999). 

We must emphasise that our models only account for radio emission produced from 
star-formation processes. The sources labelled as elliptical (or early type) in Fig. [K] are not 
expected to lie on any of the star-formation derived locii. An AGN most likely dominates 
their radio emission. We include them here merely as a comparision, and their relationship 
to normal starbursts is discussed further below. 

The sources in Fig. ^ appear more-or-less consistent with the dusty "0.1 Gyr 
starburst" model. This could in principle apply to the two sources with spiral/disk-like 
morphology (labelled "S"), but is unconventional for the five elliptical morphologies 
(see below). A comparison with radio-to-near-infrared flux ratios further constrains the 



underlying properties of these sources. Figure shows the radio-to-near-infrared flux 
ratio R(1A/H), defined analagous to eq.||, as a function of r — H color. The near-infrared 
emission is dominated by old stars and is less affected by dust than the optical. The 
radio-to-near-infrared flux ratio should therefore be relatively insensitive to dust. Given 
the simplicity of our models, the two disk-like sources may not necessarily represent "0.1 
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Fig. 10. — Radio-to-optical (r-band) flux ratio (see eq.[ H) as a function of r — H color 
for radio sources with available optical and near-infrared detections. The predictions of four 
synthesis models for < z < 1.5 (z — 0-starting bottom left of curves to z = 1.5-top right) 
are also shown. These assume no extinction (thin curves) and with a rest-frame extinction 
Ay = 2.5 mag (thick curves). 
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Gyr starbursts" as indicated in Fig. [I]]. They could also belong to the Sab/Sbc or Scd/Sdm 
classes. For this to be true, the following additional components may play an important 
role in more 'realistic' models: 'optically-thick' dust that completely obscures both the H 
and r emission without causing appreciable reddening in r — H color, or, contamination by 
at least an order of magnitude times more radio emission from a central AGN than that 
produced purely by supernovae. The second possibility is favored by radio observations of a 
number of luminous infrared galaxies by Norris et al. (1988), where some showed evidence 
for significant radio emission from compact Seyfert-like nuclei. 

The large discrepancy between the four sources with elliptical morphologies (labelled 
"E") and predictions from the early-type E/SO models suggests the importance of a 
significant AGN contribution to the radio emission. Appreciable amounts of optically-thick 
dust suppressing the optical and near-infrared light (except for extinction by diffuse, 
optically-thin dust) is not favored by observations (eg. Goudfrooij & de Jong 1995). Most, if 
not all of these ellipticals are likely to be radio-powered by AGN. At the limiting sensitivity 
of our radio survey (~ 0.3 mJy), a nominal FR-I galaxy (with L1.4 = 10 24 WHz _1 ) could 



be detected to z ~ 1.3, and indeed the spread in r — H color for the ellipticals in Fig. |TT 
is consistent with the E/SO (optically-thin dust) model to this redshift. A comparison 
between the model R(1A/H) values with actual observed values implies that such AGN will 
contribute a factor 10 2 times more radio emission than that produced by any underlying star 
formation activity in these systems. It is important to note that the ratio of AGN-to-stellar 
powered radio activity has a huge spread for the elliptical population in general, and that 
the factor 10 2 only illustrates a property specific to the ellipticals in our radio sample. 

To summarise, our use of a simple synthesis model that includes radio emission and 
dust reddening to analyse the properties of sub-mJy radio sources has shown the following: 
first, the presence of dust with extinctions Ay ~ 2 mag and possibly greater, consistent 
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Fig. 11. — Radio-to-near-infrared (if -band) flux ratio as a function of r — H color and 
models which include no extinction (thin curves) and A v = 2.5 mag (thick curves). 
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with previous more direct determinations, and second, that the level of radio emission from 
non-stellar processes such as AGN could be easily inferred and constrained. This will be 
particularly important for starbursts hosting Seyfert nuclei where a comparison with more 
sophisticated dust models may be required to infer the relative contributions. 



6. A Method to Select "ULIGS" via Radio/Optical Color 

Since the emission (and dust absorption) properties from normal galaxy populations 
are reasonably well known, a color-color diagram such as Fig. [H] could provide a potential 
diagnostic for selecting ultraluminous infrared galaxies (ULIGS) to high redshift. The 
relatively low sensitivity of the Infrared-Astronomical Satellite (IRAS) has primarily 
confined ULIG selection to the local Universe (Sanders & Mirabel 1996), although there 
is some speculation that recently discovered faint "SCUBA" sources at sub-millimetre 
wavelengths could represent their high-redshift counterparts (eg. Blain et al. 1999). 
Approximately 80% of local ULIGS are believed to be powered by starbursts and the 
remainder show evidence for an AGN contribution (Genzel et al. 1998; Lutz et al. 1998). 
Far-infrared observations have shown that dust and molecular gas in local ULIGS is 
concentrated in compact regions <^ lkpc (Okumura et al. 1991, Bryant 1996) and that a 
large fraction of the optical/UV emission is hidden by optically-thick dust (Sanders et al. 
1988). A study of their properties and importance to galaxy evolution therefore requires 
observations at wavelengths virtually immune to dust absorption. Radio frequencies provide 
an excellent window of opportunity. 



Figure [12| illustrates the predicted locus in color-color space using the synthetic SEDs 
of three local far-IR selected systems: Arp 220 (L IR ~ 1.6 x 10 12 Lq) - a ULIG undergoing 
a powerful starburst as seen via high resolution radio observations by Smith et al. (1998); 
M82 (L /K ~ 6 x 1O 1O L ) - a system under going a weak-to- moderate starburst, and Mrk 
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273 (Ljr ~ 2.6 x 10 12 Lq) - a ULIG whose bolometric emission is believed to be dominated 
by a hidden central AGN from the presence of strong Seyfert-2 lines and moderately strong 
hard X-ray (2-10 keV) emission (Turner et al. 1997). We have used the synthetic SEDs 
generated by Devriendt et al. (1999) to model the starburst emission. These authors used 
a self-consistent modelling approach to predict the stellar optical/UV/near-IR emission, its 
reprocessing into the mid-IR-to-sub-mm by dust, and the nonthermal stellar-powered radio 
emission based on the empirical radio-to-far-IR luminosity correlation. Due to its strong 
AGN-dominated nature, the starburst synthetic SED predicted by Devriendt et al. for Mrk 
273 differs appreciably from that observed in the radio. For this source, we therefore used 
the Devriendt et al. SED at wavelengths A < 1mm and extrapolated into the radio using 
its actual observed radio-to-lmm spectral slope and fluxes (obtained from the NASA/IPAC 
Extragalactic Database []). 



Figure [12] shows that a galaxy characteristic of the (low-IR luminous) M82 system 
will occupy a region similar to that occupied by normal galaxies in this study (and 
also their predicted synthetic colors in Fig. [10D. Luminous systems classified as ULIGS 
however (Arp 220 and Mrk 273), will tend have higher radio-to-optical flux ratios which 
could be easily selected. This can be explained by the well-observed correlation between 
far-infrared luminosity and far-IR (60 — 100/mi)-to-optical spectral slope (Soifer et al. 
1987). Consequently, the most IR-luminous systems with the largest far-IR-to-optical 
ratios are also likely to have a high level of radio-emission due to its strong correlation with 
IR luminosity. This will lead to a larger than average radio-to-optical flux ratio for ULIGS 



in general as shown in Fig. 12. 



7 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. 



-40- 



10 



8 



DC 







z=6 






o 




****** i 


o 












o 






o 






o This study 




o 


Arp 220 (ULIRG; Starburst) 






M 82 (Low L IR Starburst) 






Mrk 273 (ULIRG; AGN) 

i 




3 


5 7 



r - H (mag) 



Fig. 12. — Locus of observed radio-to-optical flux ratio versus r — H color for < z < 6 
(Az = 0.2) using three local far-IR selected systems: Arp 220 (L IR ~ 1.6 x 10 12 Lq; 
strong starburst), M82 (L IR ~ 6 x 10 10 Lq; weak/moderate starburst) and Mrk 273 
(L IR ~ 2.6 x 1O 12 L ; AGN dominated). 
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The existence of systems with either larger rest-frame optical/UV extinction or excess 
AGN contribution to the radio than the ULIGS considered here will be shifted further 
upwards on this plot. Diagnostics to distinguish between AGN and starburst dominated 
ULIGS using radio-to-optical color alone will not be trivial and is left to a future study. 
The three ULIGS in Fig. [12| represent a range of known ULIGS and their location on this 
plot simply serves as a diagnostic to pre-select ULIG candidates for further study. 

A system like Arp 220 (with uL u (lAGB.z) ~ 2.5 x 10 6 Lq) could be observed to 
redshift z ~ 1.60 if initially selected from a radio survey limited to <Si.4GHz — 50// Jy. Thus, 
to limiting sensitivities reached by existing 1.4 GHz surveys, such a method may not 
probe the highest redshifts. Nonetheless, as shown in Fig. such systems could still be 
well separated from normal galaxies to this redshift. Assuming an Arp 220-like SED and 
moderate luminosity evolution (L 60Atm oc (1 + -z) 2 ' 5 ), the surface density of ULIGS to z ~ 1.6 
is expected to be of order 150 deg -2 at Si.4gh z ^ 50/iJy, or about 6% of the integral count 
to this sensitivity (Richards 2000). They should therefore exist in significant numbers in 
deep large-area radio surveys. 



7. SUMMARY AND CONCLUSIONS 

We have used the VLA radio telescope to image a contiguous 33 x 33 arcmin 2 area to a 
(mean) limiting (5<r) sensitivity of ~ 0.35 mJy From a total of 62 detections, the results of 
optical and near-infrared photometry are reported for 43 sources. Our optical photometry 
is more sensitive than previous optical follow-up studies of radio surveys of similar depth. 
Our main findings are: 

(1) We have used a robust, likelihood-ratio method for determining optical 
8 Assumes q = 0.5, H = SOkms^ 1 Mpc -1 
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identifications and their reliability. This method is seldom used in identification studies 
and is insensitive to assumptions concerning fluctuations in background source density 
and Gaussian error distributions. We assigned optical candidates to 26 radio sources with 
reliability ^ 80%. Nine radio sources are uncertain and/or ambiguous, and eight are empty 
fields. Near-infrared photometry from the 2MASS database was reported for 7 sources. 

(2) The eight optical empty field sources all display compact and symmetric radio 
morphologies and most probably represent compact starbursts at z <^ 0.3 strongly obscured 
by dust. They may require at least 4 magnitudes optical extinction to account for their 
large radio-to-optical flux ratio compared to the identified population. Our conclusion for 
them being 'compact starbursts' is very tentative as it is purely based on starburst versus 
AGN number statistics expected from sub-mJy radio surveys. Further deep infrared/optical 
imaging and spectroscopy will be necessary. 

(3) Consistent with previous studies, our deep (r ~ 25) optical imaging shows that 
the optical appearence can be divided into two classes according to radio flux-density: 
elliptical-like morphologies for ^ 2 mJy, and peculiar or disturbed for <^ 2 mJy. 

(4) Using a stellar synthesis model which includes radio emission and dust reddening, 
we find that the near-infrared-to-optical emission in a small, bright sub-sample is reddened 
by 'optically thin' dust with A v ~ 2 — 2.5 mag, regardless of morphological type. This 
appears consistent with other more direct determinations. Consistent with previous studies, 
the radio emission from early-type systems seems to be powered by AGN rather than 
star-formation to account for their anomalously large radio-to-optical(-near-infrared) 
ratios. 

(5) Our analysis shows that a radio/optical or radio/near-IR color selection technique 
could provide a potential means for detecting ULIG-type objects to z ~ 1.6. 
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Despite the lack of spectroscopic information, our study of a homogeneous population 
of faint radio sources has stressed the importance of dust on studies of intrinsic galaxy 
properties and their evolution at optical wavelengths. A future goal would be to obtain 
spectra, or multi-color optical/near-infrared photometry to better explore these sources 
and the validity of the simple stellar synthesis models presented in this paper. The ever 
improving resolution (and sensitivity) capabilities of optical/near-IR detectors over those 
feasible at (the longest) radio wavelengths requires robust identification techniques to better 
ascertain their properties. Likelihood-ratios provide one such technique. The present study 
complements other deep optical studies of faint radio sources to constitute a statistically 
significant sample for inferring their nature and importance to galaxy evolution. 
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Table 1. The regions used for the source extractions 



Region 


RA(J2000) a 


Dec(J2000) a 




5cr limit 


# sources b 








(mJy) 


(mJy) 




1 


00 ft 13 m 48 s .43 


+25°48'07 /, .7 


0.061 


0.305 


15 


2 


00 ft 12 m 34 s .09 


+25°47 , 11".7 


0.065 


0.325 


15 


3 


00 /l 12 m 33Ml 


+26°03'27".7 


0.090 


0.450 


20 


4 


00 fe 13 m 48 s .21 


+26°02'47".7 


0.076 


0.380 


14 c 



a Defines the center of each region with 16.5 x 16.5arcmin . 

b Total number of sources with flux density > 5a rms . 

c Includes the separate components of two double-component sources. 



TABLE 2 
Catalog of Radio Source Data 





RA(J2000) 
(h,m,s) 


Dec 

1 
\ 


(J2000) 
"") 


<*RA 
(") 


(") 


S P 

Sm Tvl 

\ LLL ^y } 


(mJy) 




(mJy) 


#M 
(") 


(") 


PA 
(dee) 


"e M 
(") 


»8 m 

(") 


idee) 


RADJ001158+2550 


00 


[ 1 


1st 


,0o8 


+25 


oO 


24 


.41 





.37 


0. 


43 


1.468 





.095 


1.410 





.088 


15 


944 


1 3 


.579 


169 


1 


.02 / 





J 


16 


RADJ001158+2556 


00 


1 1 


58 


.078 


+25 


06 


48 


.20 





24 


0. 


26 


2.176 





,09o 


2.820 





075 


14 


482 


1 2 


.969 


24 





6.' iO 





o64 


16 


RADJ001203+2544 


00 


1 2 


03 




+25 


II 


5-1 


.18 


1 


.35 


1. 


06 


0.424 





. 084 


0.284 





.116 


15 


249 


9. 


908 


63 


3 


.402 


2 


.210 


20 


RADJ001204+2540 




'"' 


n! 


fin 


+25 






21 








Vt 


2.653 




nil* 


2.930 




' ' 


If 


™ 


| 3 




51 


° 


^' 




^ 


9 


RADJ001204+2605 


00 


1 2 


04 


.572 


+26 


Or? 


1 


.42 





66 


0, 


6 1 


0.973 





.084 


0.962 





0. i6 


18 


466 


1 2 


.055 


129 


1 


.795 


1 


172 


9 


RADJ001206+2611 


00 


1 2 


06 


.178 


+26 


1 1 


05 


.20 





.49 


0. 


50 


1.096 





.095 


1.050 





.046 


16 


4,30 


1 3 


.115 


19 


1 


.419 


1 


.132 


15 


RADJ001209+2540 


00 


1 2 


Of) 


.614 


+25 


40 


24 


.41 





.30 


0. 


3 1 


1.880 





.09,3 


2.780 





.099 


14 


993 


1 4 


.237 


29 





.755 





.717 


39 


RADJ001210+2551 


00 


1 2 


10 


.38 1 


+25 


5 1 


36 


.54 





.40 


0. 


43 


1.402 





.094 


1.300 





110 


1 6 


091 


1 2 


,989 


147 


1 


086 





876 


13 


RADJ001216+2601 


00 


1 2 


1 


.041 


+26 


1 


.3 1 


.64 





03 


0. 


03 


17.676 





.064 


18.42 





.108 


17 


051 


1 5 


.245 


10 





.091 





.081 


2 


RADJ001216+2543 


00 


1 2 


1 6 


.854 


+25 


43 


02 


09 





60 


0. 


57 


0.955 





.084 


0.836 





140 


1 4 


521 


1 .3 


.550 


100 


1 


437 


1 


341 


58 


RADJ001217+2602 


00 


1 2 


1 7 


.177 


+26 


02 


53 


.06 





73 


0. 


6 1 


0.878 





.044 


0.844 





.159 


16 


079 


1 3 


.462 


99 


1 


.732 


1 


450 


24 


RADJ001218+2554 


00 


1 2 


1 8 


.056 


+25 


54 


44 


.09 





.13 


0. 


13 


4.580 





.095 


4.700 





.173 


15 


,30 1 


1 5 


.08 1 


47 





. 3 1 6 





,3 1 1 


18 


RADJ001218+2539 


00 


1 2 


1 8 


.193 


+25 


30 


42 


.82 





.07 


0. 


07 


8.443 





.095 


8.380 





.167 


14 


927 


1 4 


.783 


16 





.167 





166 


17 


RADJ001227+2543 


00 


1 2 


27 


.184 


+25 


■ 1 3 


30 


. 1 


3 


.09 


3. 


67 


0.240 





. 090 


0.408 





207 


26 


143 


1 4 


.62.3 


144 


9 


884 


5 


529 


25 


RADJ001230+2610 










+26 
















0.734 






1.141 














134 










5 


RADJ001237+2607 


00 


1 2 


37 


254 


+26 


07 


48 


87 


o 


54 




53 


1.157 


o 


089 


1.230 


o 


147 


16 


8^3 


1 ^ 


221 


49 


1 


375 


^ 


163 


20 


RADJ001240+2609 


00 


12 


40 


.394 


+26 


09 


57 


.71 


1 


.18 


1. 


57 


0.447 





.094 


0.419 





.157 


18 


876 


1 1 


.192 


153 


3 


,994 


2 


,368 


16 


RADJ001241+2610 


00 


12 


1 1 


.007 


+26 


10 


.31 


.16 





.31 


0. 


28 


2.006 





.071 


1.970 





. 1 16 


16 


611 


13 


.268 


121 





.783 





.626 


8 


RADJ001242+2600 


00 


12 


42 


.209 


+26 


00 


33 


.55 





.18 


0. 


18 


3.229 





.095 


4.160 





.086 


15 


278 


14 


.392 


62 





.448 





.422 


20 


RADJ001243+2558 


00 


12 


43 


.323 


+25 


58 


10 


.05 





.51 


0. 


50 


1.158 





.095 


1.130 





136 


15 


394 


14 


.207 


52 


1 


.257 


1 


160 


41 


RADJ001246+2543 


00 


12 


46 


.030 


+25 


43 


51 


.08 


2 


.24 


1. 


63 


0.344 





093 


0.410 





123 


19 


944 


13 


.441 


107 




.430 


3 


659 


17 


RADJ001247+2602 


00 


12 


47 


.680 


+26 


02 


38 


.49 





.86 


0. 


00 


0.559 





.094 


0.377 





.127 


12 


744 


1 1 


.914 


21 


2 


.155 


2 


.014 


12 


RADJ001250+2550 


00 


12 


.",0 


.036 


+25 


50 


10 


.68 





.45 


0. 


46 


1.224 





.095 


1.050 





.151 


14 


153 


1.3 


.629 


26 


1 


.094 


1 


.054 


83 


RADJ001252+2549 


00 


12 


52 


.006 


+25 


40 


38 


.49 


1 


.21 


0. 


88 


0.655 





.033 


0.855 





.130 


20 


303 


14 


.466 


82 


2 


.881 


2 


.053 


17 


RADJ001254+2600 


00 


12 


54 


.835 


+26 


00 


48 


.52 





.34 


0. 


34 


1.656 





.095 


1.490 





.125 


14 


504 


13 


.978 


130 





.829 





.799 


63 


RADJ001255+2541 


00 


12 


55 


.365 


+25 


41 


53 


.48 





14 


0. 


12 


4.717 





.094 


5.470 





.182 


17 


566 


14 


.852 


96 





351 





.296 


5 


RADJ001257+2559 


00 


12 


57 


.288 


+25 


50 


43 


.59 





.84 


1. 


33 


0.544 





.093 


0.553 





.155 


18 


083 


1 1 


.514 


179 


3 


.141 


2 


.000 


15 


RADJ001259+2607 


00 


12 


59 


.107 


+26 


07 


18 


.91 


1 


.86 


1. 


44 


0.497 





.089 


0.501 





.218 


24 


476 


18 


.828 


84 


4 


.406 


3 


.390 


27 


RADJ001300+2556 


00 


13 


00 


.983 


+25 


56 


57 


.10 





18 


0. 


18 


3.230 





.095 


4.070 





.096 


15 


100 


14 


.147 


129 





.442 





.414 


18 


RADJ001303+2602 


00 


13 


03 


.185 


+26 


02 


07 


.41 





11 


0. 


1 1 


5.136 





.095 


5.270 





.097 


15 


289 


15 


.108 


112 





.282 





,278 


23 


RADJ001304+2544 


00 


13 


04 


.017 


+25 


44 


25 


.38 





.89 


0. 


96 


0.545 





.094 


0.977 





.129 


1.3 


924 


1 1 


.175 


146 


2 


.416 


1 


,939 


32 


RADJ001305+2607 


00 


13 


05 


.397 


+26 


07 


54 


.15 





.18 


0. 


16 


3.906 





.093 


4.990 





.119 


19 


377 


14 


.842 


115 





.462 





,354 


4 


RADJ001307+2542 


00 


13 


07 


.094 


+25 


42 


14 


.25 





.25 


0. 


25 


2.384 





.095 


3.380 





.091 


15 


068 


14 


.900 


7 





.598 





,591 


13 


RADJ001311+2556 


00 


13 


1 1 


.181 


+25 


56 


1 1 


.98 





.12 


0. 


12 


4.782 





.095 


4.900 





.076 


14 


915 


14 


.203 


150 





.295 





.281 


16 


RADJ001313+2545 


00 


13 


13 


.019 


+25 


45 


15 


.15 


1 


.57 


1. 


75 


0.358 





.094 


0.529 





155 


17 


375 


1 1 


.916 


36 


4 


.588 


3 


.147 


28 


RADJ001314+2541 


00 


13 


M 


.646 


+25 


41 


M 


.85 





.74 


0. 


80 


0.723 





.074 


0.824 





121 


14 


606 


1.3 


.296 


162 


1 


.909 


1 


.738 


16 


RADJ001318+2604C1 


00 


13 


IS 


.834 


+26 


04 


38 


.76 


1 


.28 


1. 


09 


0.578 





.093 


0.597 





.156 


22 


405 


9. 


339 


128 


3 


666 


1 


.528 


7 


RADJ001320+2604C2 


00 


13 


20 


.122 


+26 


04 


1 1 


.43 





.05 


0. 


05 


12.512 





.092 


13.44 





.081 


17 


148 


14 


.930 


130 





.129 





.113 


2 


RADJ001320+2556 


00 


13 


20 


.552 


+25 


56 


55 


.59 





.58 


0. 


71 


0.911 





.094 


0.993 





154 


17 


.319 


1 1 


.884 


151 


1 


.797 


1 


.233 


1 1 



1 



TABLE 2— Continued 



Name RA(J2000) Dec(J2000) <j RA <JDec S P rr Sp S T rrs T 6m 9m PA <Je M <re„ "PA 

(h,m, S ) ("'") (") (") (mJy) (mjy) (mjy) (mjy) (") (") (deg) (") (") (deg) 



CD 



RADJ001321+2543 


00 


13 


21 


.405 


+25 


43 


49 


.37 





.69 


0. 


70 


0.800 





.096 





.888 





148 


14 


704 


13 


.115 


140 


1 


.745 


1 


551 


10 


RADJ001323+2600 


00 


13 


23 


.522 


+26 


00 


29 


.25 





.03 


0. 


03 


15 


..877 





.095 


16.24 





.097 


15 


258 


15 


.080 


130 





.091 





.090 


21 


JrvAJJ JUUloiO+iiooo 


























1. 


170 






1 


.700 


























RADJ001327+2543 


00 


13 


27 


.289 


+25 


43 


35 


.07 


1 


.29 


1. 


22 


0. 


424 





.094 





.332 





119 


13 


870 


12 


.706 


64 


3 


.091 


2 


.832 


13 


RADJ001339+2543 


00 


13 


39 


.899 


+25 


43 


54 


.55 





18 


0. 


18 


3. 


140 





.065 


4 


.010 





161 


14 


983 


14 


.384 


141 





.451 





.433 


30 


RADJ001342+2546 


00 


13 


42 


.203 


+25 


46 


IS 


.76 


2 


.95 


1. 


55 


0. 


361 





.030 





.654 





.237 


27 


823 


II 


.628 


92 


6 


.953 


3 


.656 


15 


RADJ001349+2539 


00 


13 


49 


.176 


+25 


39 


20 


.08 





.42 


0. 


35 


1. 


519 





.035 


1 


.450 





.162 


16 


018 


13 


.413 


88 





.998 





.835 


14 


RADJ001349+2609 


00 


13 


49 


.935 


+26 


09 


25 


.19 





10 


0. 


10 


6. 


301 





.094 


7 


.330 





.088 


17 


944 


14 


.591 


46 





.268 





.218 


3 


RADJ001350+2608 


00 


13 


50 


.742 


+20 


08 


06 


.93 


1 


.31 


1. 


04 


0. 


460 





.094 





.364 





.110 


15 


341 


1 1 


.627 


107 


3 


.153 


2 


.390 


30 


RADJ001355+2540 


00 


13 


55 


.638 


+25 


40 


48 


.82 


1 


.12 


1. 


15 


0. 


538 





.094 





.547 





165 


16 


771 


13 


.643 


138 


2 


.945 


2 


.396 


34 


RADJ001355+2553 


00 


13 


55 


.874 


+25 


53 


27 


.52 





.73 


0. 


72 


0. 


760 





.097 


1 


.432 





135 


15 


362 


12 


.182 


133 


1 


.912 


1 


.516 


22 


RADJ001359+2601 


00 


13 


59 


.764 


+26 


01 


23 


.68 





.41 


0. 


39 


1. 


413 





.095 


1 


.360 





.045 


15 


092 


13 


.286 


124 


1 


.01 1 





.890 


21 


RADJ001359+2605 


00 


13 


59 


.945 


+26 


05 


17 


.72 





.49 


0. 


63 


1. 


080 





.095 


1 


.090 





.087 


17 


033 


13 


.351 


175 


1 


.492 


1 


.169 


14 


RADJ001401+2559C1 


00 


14 


01 


.544 


+25 


59 


59 


.45 





.06 


0. 


07 


9. 


410 





.094 


9 


.860 





.132 


17 


264 


15 


.048 


150 





173 





151 


3 


RADJ001401+2610 


00 




01 


.825 


+26 


10 


59 


.96 


1 


.14 


1. 


04 


0. 


465 





.084 





.329 





.120 


13 


499 


1 1 


.792 


114 


2 


.742 


2 


.395 


31 


RADJ001403+2559C2 


00 




03 


.162 


+25 


59 


37 


.11 





08 


0. 


08 


8. 


005 





.088 


9 


080 





.028 


16 


501 


15 


.463 


140 





195 





.182 


7 


RADJ001405+2551 


00 


11 


05 


.675 


+25 


51 


29 


.74 





37 


0. 


30 


1. 


884 





.095 


2 


100 





.081 


17 


524 


14 


.319 


90 





.879 





.719 


10 


RADJ001408+2550 


00 


14 


08 


.057 


+25 


50 


37 


.67 





.50 


0. 


47 


L. 


172 





.035 


1 


.650 





.072 


14 


896 


13 


.522 


113 


1 


.202 


1 


.091 


24 


RADJ001410+2610 


00 


14 


10 


.900 


+26 


10 


35 


.97 





.32 


0. 


30 


1. 


828 





.095 


1 


.870 





.094 


15 


161 


13 


.586 


60 





.785 





.703 


19 


RADJ001411+2546 


00 


14 


1 1 


.017 


+25 


46 


03 


.46 





.23 


0. 


20 


2. 


919 





.095 


3 


.190 





1 10 


16 


703 


II 


.714 


93 





.541 





.477 


10 


RADJ001411+2541 


00 


14 


1 1 


.850 


+25 


4 1 


24 


.11 


1 


.19 


1. 


40 


0. 


417 





.088 





.331 





.140 


14 


SI 1 


12 


.073 


19 


3 


.356 


2 


.736 


45 


RADJ001416+2610 


00 


14 


16 


.812 


+26 


10 


51 


.50 





.86 


0. 


94 


0. 


681 





.093 





.696 





.126 


16 


754 


13 


.742 


149 


2 


.327 


1 


.909 


28 


RADJ001418+2540 


00 


14 


18 


.769 


+25 


40 


08 


.27 





.49 


0. 


54 


1. 


107 





.095 


1 


.200 





.110 


15 


381 


13 


.182 


150 


1 


.314 


1 


.126 


12 



2 



TABLE 3 
Optical-Infrared Data 



Name 


RA(J2000) a Dec(J2000) a S rad - op i \ag(LR) R, d s 
(h,m ]S ) (*'") (") 


r 




H(o H ) 


K(o K ) 


r — H r — K morph c 















RADJ001216 + 2601 d 


00 


12 


in 


530 


+26 


m 


32 


77 


7 


41 


-4.163 


< 0.1 


24.0S3(? 


RADJ001216 + 2543 


00 


12 


in 


783 


+25 


43 


m 


34 


1 


31 


3.537 


=3 1 


17.565 


RADJ001217+2602 


00 


12 


17 


280 


+26 


02 


54 


85 


1 


78 


1.514 


=3 1 


22.057 


RADJ001218+2554 


























> 25 


RADJ001227+2543 


00 


12 


27 


033 


+25 


43 


41 


19 


3 


09 


0.685 


0.862 


22.971 


RADJ001237+2607 


00 


12 


37 


248 


+26 


07 


48 


49 





39 


2.883 


sa 1 


19.371 


RADJ001242 + 2600 d 


00 


12 


42 


347 


+26 


00 


30 


S3 


3 


42 


0.278 


0.783 


23.621(? 


RADJ001243+255S 


























> 25 


RADJ001 246 + 2543 


00 


12 


45 


777 


+25 


43 


50 


59 


3 


83 


1.293 


S3 1 


21.169 


RADJ001247+2602 d 


00 


12 


47 


880 


+26 


02 


35 


96 


4 


03 


0.302 


0.788 


22.76S(? 


RADJ001250 + 2550 


00 


12 


50 


137 


+25 


50 


10 


86 


1 


52 


1.225 


0.995 


23.128 


RADJ001252 + 2549 


00 


12 


r.i 


9 f'i 


+25 


49 


38 


81 





65 


2.919 


=3 1 


18.984 


RADJ001254+2600 d 


00 


12 


54 


855 


+26 


00 


44 


70 


3 


S3 


0.246 


0.764 


23.004(? 


RADJ001255 + 2541 d 


00 


12 


55 


064 


+25 


■1 1 


54 


13 


4 


56 


-0.689 


< 0.1 


23.599(? 


RADJ001257+2559 


00 


12 


57 


287 


+25 


59 


■13 


60 





03 


2.248 


=3 1 


20.384 


RADJ001259 + 2607 


00 


12 


59 


187 


+'2 f'i 


07 


19 


09 


1 


21 


2.706 


=3 1 


19.010 


RADJ001300 + 2556 


00 


13 


01 


007 


+25 


56 


57 


08 





35 


1.505 


=3 1 


22.974 


RADJ001303+2602 


00 


13 


03 


204 


+2 f'i 


02 


07 


27 





31 


1.681 


=3 1 


22.492 


RADJ001304+2544 


























> 25 


RADJ001305 + 2607 


00 


13 


05 


293 


+26 


07 


53 


81 


1 


60 


1.675 


=3 1 


21.835 


RADJ001307+2542 


























> 25 


RADJ001311 + 2556 


00 


13 


1 1 


fifi 3 


+25 


56 


41 


19 


1 


95 


1.515 


=3 1 


21.950 


RADJ001313+2545 


00 


13 


13 


049 


+25 


45 


17.44 


2 


33 


0.778 


0.882 


23.464 


RADJ001314+2541 


00 


13 


14 


595 


+25 


41 


14 


73 





78 


1.118 


0.967 


23.698 


RADJ00131S+2604C1 


00 


13 


1 8 


794 


+2 f'i 


04 


37 


24 


1 


64 


1.189 


0.984 


22.847 


RADJ001320+2604C2 


00 


13 


20 


IfsO 


+2 f'i 


04 


10 


69 


1 


37 


1.706 


=3 1 


21.938 


RADJ001320 + 2556 


00 


13 


20 


705 


+25 


56 


53 


85 


2 


88 


1.647 




20.481 


RADJ001321 + 2543 


00 


13 


21 


339 


+25 


43 


48 


S 1 


1 


15 


0.996 


0.936 


22.929 


RADJ001323+2600 d 


00 


13 


23 


256 


+26 


00 


27 


73 


4 


28 


0.247 


0.775 


24.040(? 


RADJ001325 + 2553 


00 


13 


25 


803 


+25 


53 


39 


14 


2 


75 


0.695 


0.864 


23.650 


RADJ001327+2543 


00 


13 


27 


260 


+25 


43 


33 


55 


1 


58 


1.265 


=3 1 


22.639 


RADJ001339 + 2543 


























> 25 


RADJ001342 + 2546 d 


00 


13 


41 


666 


+25 


46 


in 


96 


8 


27 


0.238 


0.770 


21.194(? 


RADJ001350 + 260S 


00 


13 


50 


572 


+2 f'i 


08 


12 


05 


4 


73 


0.617 


0.858 


22.178 


RADJ001355 + 2553 


00 


13 


55 


857 


+25 


53 


28 


22 





75 


3.729 


=3 1 


17.286 


RADJ001359 + 2601 


























> 25 


RADJ001359 + 2605 d 


00 


13 


59 


795 


+26 


05 


22 


01 


4 


84 


-0.760 


< 0.1 


23.763(? 


RADJ001401+2559C1 


00 


1-1 


n i 


559 


+25 


59 


57 


15 


2 


3 1 


0.861 


0.903 


23.390 


RADJ001403+2559C2 d 


00 


14 


03 


OSS 


+25 


59 


33 


06 


4 


20 


-0.790 


< 0.1 


24.707(? 


RADJ001405 + 2551 


00 


14 


05 


728 


+25 


r.i 


28 


88 


1 


16 


1.893 


=3 1 


21.508 


RADJ00140S+2550 


























> 25 


RADJ001411 + 2546 


























> 25 


RADJ001411 + 2541 


00 


14 


1 1 


8 05 


+25 


41 


22 


34 


1 


93 


1.016 


0.938 


23.164 



I H. 05(1(0.003) 15.409(0.138) I 4.525 (0. 08ii) 



17.416(0.285) 16.057(0.204) 



17.066(0.210) 
15.158(0.074) 



16.131(0.225) 
14.290(0.075) 



16.118(0.207) 



> 15.178 
13.913(0.074) 



.128) 15.224(0.179) 7 



15.642(0.071) 14.756(0.082) 14.419(0.098) 



927 <2.f"i63 



253 <5 
720 



206 
097 



393 7.7 



'"'TVf.TR to optical position. 

b Values with ">" are r la upper limits on measured fluxes, otherwise, errors are quoted in parenthesis. 
c Optical morphology: E — Elliptical, Irr — Irregular or Disturbed, S — Spiral, ? — Unknown. 

d Optical identification uncertain due to either extended radio emission and large offset from radio centroid position, or, 'moderately low' LR value. 
e Reliability of the optical identification. Cases with 'very low' and 'very Inrgo loqiT.Ti.) mv listed with Rid ^ 1 and Rid ^ 1 respectively. 
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Table 4. Population synthesis models 



Type 


SFR(t) 


IMF 


Age (Gyr) a 


E/SO 


rf 1 exp (-t/n) 


Scalo 


12.7 


Sab/Sbc 


r 8 _1 exp (-t/r 8 ) 


Scalo 


12.7 


Scd/Sdm 


constant 


Salpeter 


12.7 


SB 


constant 


Salpeter 


0.1 



a For the SB-type, an age of 0.1 Gyr is assumed to 
apply at all redshifts. For other types, this refers to 
the present day age. 
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